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SUMMARY 


A method is presented which allows the use of nonlinear section 
lift data in the calculation of the spanwise lift distribution of 
unswept wings with flaps or ailerons. This method is based upon lifting- 
line theory and is an extension to the method described in NACA Rep. 865. 
The mathematical treatment of the discontinuity in absolute angle of 
attack at the end of the flap or aileron involves the use of a correction 
factor which accounts for the inability of a limited trigonometric series 
to represent adequately the spanwise lift distribution. ' A treatment of 
the apparent discontinuity 'in maximum section lift coefficient is also 
described. In order to minimize the computing time and to illustrate 
the procedures involved, simplified computing forms containing detailed 
examples are given for both symmetrical and asymmetrical lift distribu- 
tions. A few comparisons of calculated characteristics with those 
obtained experimentally are also presented. 


INTRODUCTION 


Unswept-wing characteristics calculated by the method of refer- 
ence 1, in which nonlinear section lift data are used, have been found 
to agree much closer with experimental data in the region of maximum 
lift coefficient than those calculated by methods in which linear section 
lift curves -are used. It appears feasible that the similar use of non- 
linear section data would yield improved results for unswept wings with 
flaps or ailerons. The deflection of a partial-span flap or aileron, 
however, causes discontinuities in the spanwise distribution of the 
absolute angle of attack. If such discontinuities exist, an excessively 
large number of spanwise stations must be considered in order to obtain 
a solution by the method of reference 1 or by any other method in which 
the lift distribution is approximated by a trigonometric series. A 
different method of treatment of the discontinuity is therefore desirable. 
Developed herein is a new method of treatment involving the use of a 
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correction factor which accounts for the inability of a limited trigo- 
nometric series to represent adequately the spanwise lift distribution 
of a wing with partial-span flaps or ailerons deflected. This correc- 
tion factor is obtained with the aid of reference 2 and is used in con- 
junction with the system of multipliers of reference 1 to obtain the 
induced angle of attack from the spanwise lift distribution. The sub- 
sequent calculation of the lift distribution by means of successive 
approximations is similar to that of reference 1. 

The mathematical treatment of the discontinuity in the spanwise 
distribution of absolute angle of attack is only part of the problem 
involved in calculations for wings with flaps or ailerons. The direct 
use of two-dimensional. lift data would indicate a discontinuity in the 
spanwise distribution of maximum section lift coefficient. Obviously, 
the flow about the wing sections near the end of a flap or aileron is 
not two-dimensional. A rational method of obtaining three-dimensional 
section data from the two-dimensional data has therefore been devised 
and is presented herein. This method is substantiated by experimental 
pressure distributions. 

In addition to the presentation of the method of making the calcu- 
lations, simplified computing forms are given and their use is illus- 
trated by detailed examples for both symmetrical and asymmetrical dis- 
tributions. A few comparisons of calculated results with experimental 
data are also given. 

This paper and reference 1 are intended to supplement each other. 
The reader is therefore expected to be reasonably familiar with 
reference 1. 


SYMBOLS 


As used herein, the term "section" designates the characteristic 
of a section in three-dimensional flow. 

A aspect ratio. 

A n coefficients of trigonometric series for lift distribution 

Cp^ induced drag coefficient 

Cl wing lift coefficient 

C 1 rolling -moment coefficient 

C n ^ induced yawing-moment coefficient 
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E 

E' 


F 

R 

S 

a o 

b 

c 

C S 

c 

c *i ' 

c l 
Ac 1 

c imax 


edge-velocity factor for symmetrical part of lift distribution 

edge-velocity factor for antisymmetrical part of lift 
distribution 

factor used in altering two-dimensional lift curves 
Reynolds number 

ssure - Local static pre 
Dynamic pressure 

section lift -curve slope per degree 
span of wing 
local chord of wing 
root chord 

mean geometric chord (b/A) 
section induced-drag coefficient 
section lift coefficient 
stall margin 

maximum section lift coefficient 


ssure\ 


pressure coefficient 


^ Total pre 


(ci ) maximum two-dimensional lift coefficient 
\ t max( ) • 


Ac- 


max 


increment in maximum lift coefficient due to flap deflection 
(two-dimensional data) 


c 1 * section lift coefficient at either end of flap or aileron 

Cj section lift coefficient for part of lift distribution 

involving no discontinuity in angle of attack 

c 1 section lift coefficient for part of lift distribution due to 

discontinuity in angle of attack 

C 7 section lift coefficient for additional lift distribution 

al 

c n section normal -force coefficient 



NACA TN 2283 


wing-tip helix angle generated by rolling wing 
intervals used in integration 
airfoil thickness -chord ratio 
spanwise coordinate 

spanwise coordinate at either end of flap or aileron 
angle of attack, degrees 

correction for induced angle of attack, degrees 
effective angle of attack, degrees 
correction for effective angle of attack, degrees 
induced angle of attack, degrees 

induced angle of attack for part of lift distribution 
involving no discontinuity in angle of attack, degrees 

induced angle of attack containing discontinuity, degrees 

angle of attack for zero lift, degrees 

angle of attack for two-dimensional lift curves, degrees 

angle of attack of root section, degrees 

uncorrected induced angle of attack, degrees 

multiplier for induced angle of attack for asymmetrical 
distributions 

multiplier for induced angle of attack for antisymmetrical 
di s tr ibut ions 

magnitude of discontinuity, in absolute and induced angles of 
attack, degrees 

angle of twist, negative if washout, degrees 
angle of twist at wing tip, degrees 
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e average angle of twist 

f '5 faired angle of twist due to flap deflection 

“5 average angle of twist due to flap deflection 

T) ratio of actual two-dimensional lift- curve slope to 

theoretical value of rt ^/90 

T] m area multiplier for asymmetrical distributions 

area multiplier for symmetrical distributions 

Q = cos - "** 

b 

e* = cos- 1 

D 

taper ratio fe - P / h ° rd J 
\R 00 t chord/ 

•multiplier for induced angle of attack for symmetrical 
distributions 

interpolation multiplier 

moment multiplier for asymmetrical distributions 
moment multiplier for antisymmetrical distributions 
used as superscript to denote value at end of flap or aileron 

DEVELOPMENT OF METHOD 
Lift Distribution 


X 

\bk 


r m 


a m 


ma 


The method involving the use of multipliers to obtain the induced 
angle of attack from the spanwise lift distribution was presented in 
reference 1. In this development, the lift distribution was approximated 
by a finite trigonometric series 


( 2 £) 

v b - fci 


(1) 
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where cos — = and m = 1, 2, 3> • • . r - 1. From the values 
r b 7 7 . 

of cjc/b for each value of m, the induced angle of attack could then 
be obtained at the points krt/r by the relation 


o-i. 



( 2 ) 


where k = 1, 2, 3> • • .r-1 and P Tn - k . denotes the multipliers for 

asymmetrical distributions. The corresponding multipliers, 

and are used for symmetrical and antisymmetrical distributions, 

respectively. These multipliers are tabulated in reference 1 for r= 20 
The method can be used directly so long r.s there is no discontinuity in 
the spanwise distribution of absolute angle of attack caused by the 
deflection of a partial-span flap or aileron. 

Lifting-line theory requires that a discontinuity in the distribu- 
tion of absolute angle of attack must be accompanied by an identical, 
discontinuity in the distribution of induced angle of attack in order to 
avoid a discontinuity in the spanwise lift distribution. An analytic 
expression for the lift distribution associated with a discontinuity in 
induced angle of attack is presented in reference 2. The complete lift 
distribution can thus be expressed as the sum of two distributions 



(3) 


where c^cybb is the distribution due to a unit discontinuity in the 
induced angle of attack and cj^cy^b is the remainder of the lift dis- 
tribution. These distributions are illustrated in figure 1. Since no 
discontinuity is associated with the distribution cj^c/b, the corre- 
sponding induced angle of attack can be obtained by means of the 
multipliers 



By definition, over the flap span 


(4) 



(5a) 
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and over the unflapped span 


C l2 = 0 


(5b) 


The total Induced angle of attack is 

a ik = a il k + a i2 k 

If, however, the multipliers were used with the total distribution 

E-l / v. r-1 /<■», r-1 

b / p 
/m 


( 6 ) 


r-x 

r 


m=l 


m 


m 


m=l ' 7 m=l 

and a 'i2^ vere added to both sides of this equation, the result would be 




Pmk 



A comparison of equation (7), for a wing with a flap or aileron, with 
equation (2), for a wing without flap or aileron, shows the addition of 
a term which is proportional to the magnitude of the discontinuity and 
acts as a correction factor to account for the inability of a limited 
trigonometric series to represent adequately the spanwise lift distribu- 
tion of a wing with partial-span flaps or ailerons deflected. For 
simplicity, equation ( 7 ) may be written as 


“Ik = °u k + 6 



( 8 ) 
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where the uncorrected induced angle of attack is expressed as 



and the correction factor per unit discontinuity is given hy the equation 



The distribution. c^c/bb given in reference 2 may be expressed in 
the form 


1 - cos(0 + 0*) rt0*sin 0 


cos (0 - 0 ) 


for = 1, 0° < 0 < 0*> for =0, 9 * < 6 < l80°; and 0* 

and 6 are in degrees. This distribution is dependent on only the 
spanwise position of the discontinuity and is independent of aspect 
ratio and taper ratio. The correction factor per unit discontinuity 
a Ck y6 given by equation (10) is therefore a function of only the span- 

wise position of the discontinuity. The distribution of ai 2 /b, for 
which equation (ll) applies, is illustrated in figure 2(a). Values of 
cj 2 c/b5 corresponding to this type of distribution are presented in 

table I for various values of 2y/b and 2y*/b. These values axe 
plotted against 2y/b in figure 3 for even increments in 2y /b and 

against 2y*/b in figure 4 for even increments in 0 = cos -1 In 

table II axe given values of a^/b corresponding to the values of 
cj 2 c /bS of table I. 

The distribution of <Xj. 2 /& illustrated in figure 2(a) would be 
applicable only for a wing with one outboard flap or aileron deflected. 
For a wing with symmetrical inboard flaps, the distribution would be 
obtained as illustrated in figure 2(b) and the values of Oc/b would 
be obtained in like manner. For example, for a wing with flaps 

extending from = 0.6 to = -0.6, the values of a^/b to be 

b b 2y* 

used would be obtained by subtracting the values of Oc/b for - 4 — =0.6 
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from those for = -0.6. Since the resulting distribution would be 

symmetrical, the values of oc/& would also be symmetrical. Similarly, 
an antisymmetrical distribution of and the corresponding values 

of clc/ 5 would be obtained as illustrated in figure 2(c). 

It should be noted that two values of ctc/5 exist at the end of 

the flap, one corresponding to + 0 (flap side of 2y* /b) and one 

2v* ° 

corresponding to - 0 (unflapped side). Only the values corre- 
sponding to + 0 are given in table II. For the unflapped side, 

b b 

These values have no practical significance unless the discontinuity 

occurs at one of the stations ^ = cos — . At these stations, either 

b r * 

of the two values may be used so long as the value is used with the 
proper section lift curve. 


In the calculation of the induced angle of attack by equation (8), 
the values of 0^/6 must be multiplied by the magnitude of the discon- 
tinuity 5. The value of 8 to be used is obtained from the section lift 
curves at the lift coefficient cj . If the discontinuity occurs at one 

of the spanwise stations ^ = cos the value of (cjc/b)* is obtained as 

one of the values of c^c/b computed. If the discontinuity occurs at 
some other position, the value of (cjc/b)* must be interpolated. This 
interpolation may be accomplished in the following manner. Even though 
the spanwise lift distribution is continuous through this point, the 
point is singular. Its singularity is due to the singularity of the 
corresponding point of the distribution of c^c/bB which has an infinite 

slope but zero radius of curvature. At this point, from equation (11), 


C *2 C f Tt0*sin 0* 
b8 / ~ BlOO 


( 12 ) 


Since both the distribution of cjc/b and the distribution of ci^cjbb con- 
tain the same type of singularity, a curve with c^c/b as the dependent 


variable and 


c 7 c 

l 2 


1 


bb as the independent variable would have no singular 
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point and could be approximated by a polynomial for which Lagrange's 
interpolation formula (reference 3) is applicable. From Lagrange's 
formula 


m'-zm.-. 


(13) 


where 


n+i 


k^m 

n+i 

_J 

k=i 

k^m 

The number (n + 1) of terms retained determines the degree (n) of the 
polynomial used in the approximation. It has been found that using 
four terms (two on each side of the desired point), which define a 
third-degree curve, gives values, at the end of the flap which agree 
within about 1 percent with the values obtained by means of the method 
of reference 4 but with about 4-5 terms of a trigonometric series. 

Values of the interpolation multipliers determined in this manner are 
presented in table III for various spanwise positions of the end of a 
flap. This type of interpolation is suitable only if the variables are 
single -valued in the range of the interpolation. This limitation pre- 
cludes the determination of multipliers for = 0.1 but this position 

is out of the range of practical flap spans. 




The method of determining the lift distribution by means of Succes- 
sive approximations is the same as that of reference 1 with the added 
step of determining the correction factor a c . For a given geometric 
angle of attack, the lift distribution is assumed. The uncorrected 
induced angle of attack is determined by equation ( 9 ) for asymmetrical 
distributions or by a corresponding equation using for symmetrical 

distributions. The multipliers 7 ^ for anti symmetrical distributions 

can be used only if the lift curves are linear and the value of 6 is 
independent of lift coefficient. The value of (c^c/b)* is obtained by 
means of equation*(l3) and is divided by (c/b)* to determine C 2 *. At 
this value of C 2 , & is read from the section lift curves as the 
difference between the effective angle of attack for the section without 
the flap and that for the section with the flap, both at the station at 
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the end of the flap. This value of 6 is multiplied by the appropriate 
values of a c /& from table II to obtain a^. The values of <Xq are 
added to the uncorrected values to obtain the values of induced 

angle of attack which are then subtracted from the geometric angle of 
attack to give the effective angle of attack at each spanvise station. 

At each value of effective angle of attack, the corresponding value 
of Cj is read from the appropriate section lift curve and multiplied 
by the value of ' c/b at that station to obtain a check value of c^c/b. 
If the check values do not agree with those originally assumed, a sec- 
ond approximation is made and the process is repeated. Further approxi- 
mations are made until one is found which is in agreement with the 
check values. 


Determination of Three-Dimensional Section Data 

If the two-dimensional lift data were used directly, a discontinu- 
ity in the spanwise distribution of maximum lift coefficient would be ' 
indicated at the end of the deflected flap or aileron. Inasmuch as this 
discontinuity does not exist in three-dimensional flow, some means of 
altering the two-dimensional data must be used to obtain what may be 
called three-dimensional section data. One method of alteration is 
described in reference 5 for use with linear section lift data/ In 
order to use this method, the lift distribution must be broken up into 
the basic lift distribution due to twist and the additional lift dis- 
tribution due to angle of attack. When nonlinear section lift data are 
used, the lift distribution cannot be broken up in this manner and the 
method of reference 5 is not applicable. For this reason, the method 
hereinafter described has* been developed. Although this development is 
not rigorous, the reasoning behind it is substantiated by experimental 
pressure distributions. 

The foundation of the method used herein is the calculation of two 
curves for which linear section lift data are used. One curve, the 
spanwise distribution of section lift coefficient for a wing with flaps 
but no other twist, is calculated by the method described herein but 
adapted for linear lift curves in a manner similar to that of refer- 
ence 1. The other curve is the additional lift-coefficient distribution 
(constant absolute angle -of -attack distribution) calculated as in refer- 
ence 1. These curves may be calculated for any convenient angle of 
attack and for any . convenient value for the discontinuity in. angle of 
attack at the end of the flap. A typical set of curves is shown in 
figure 5 . It is reasoned that the root section of a wing with partial - 
span flaps would be acting most nearly like that of a wing with full- 
span flaps so that the additional lift-coefficient distribution is mul- 
tiplied by a suitable constant to give the same value at the root. 
Likewise, it is reasoned that the tip section would be acting most 
nearly like that of a wing without flaps so that the additional 
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lift-coefficient distribution is multiplied by another constant to give 
the same value at the outermost station used in the computations (the 
0 . 9877 -semispan station when 10 points on the semispan are used). The 
differences between the .values on the curves are then divided by the 
difference between the values on the additional lift -coefficient curves 
at the end of the flap to obtain the factor F shown in figure 6 for 
several wings. In order to interpolate for the points on the additional 
lift-coefficient curves at the end of the flap, the multipliers n m , 
presented in table IV, have been determined. These multipliers were 
obtained by using 9 as the independent variable in equation (l^) 
instead of c^c/bb. It can be readily seen that the factor F is 

independent of the angle of attack and magnitude of the discontinuity 
used in the original computations. The factor F is also relatively 
independent of lift-curve slope so that it can be used near maximum lift 
where the lift curves are nonlinear. 


The maximum lift- coefficient values are altered by means of the 
factor F according to the relation 



(15) 


The values of c^ and a are then altered according to the equations 


^altered 

FIX 


unaltered 



(16) 



(17) 


The edge-velocity factor E is used in the same manner as in refer- 
ence 1. The value of E given in reference 6 is. however, probably 
more accurate than the ratio of semiperimeter to span used in refer- 
ence 1 . From reference 6 for unswept wings, 


E '\p? <18) 

The foregoing description of the alteration to the two-dimensional 
section data has been limited to wings with symmetrical partial-span 
flaps. For wings with deflected ailerons, a similar method based upon 
the same reasoning could readily be devised. 


/ 
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The' reasoning behind the method of altering the two-dimensional 
data is substantiated by the results shown in figures 7 and 8 and 
obtained from experimental pressure distributions for a wing of aspect 
ratio of 6.0, taper ratio of 0.5> and NACA 64-210 airfoil sections. 
Chordwise pressure distributions were obtained at six spanwise stations. 
Two spans of flaps, 0.49b and 0.51b, were tested so that the pressure 

distribution at ^ = 0.5 was obtained, in one case, just outboard of 

the end of the flap and, in the other case, just inboard of the end of 
the flap. As shown in figure 7, the pressure distributions at each sta- 
tion are very similar for the two flap spans, even at the station 0.5b/2. 
Furthermore, there is a gradual spanwise change from the type of loading 
associated with an airfoil section with flaps to the type associated 
with a plain airfoil. 

For comparative purposes, the section normal-force coefficients, 
obtained by integration of the chordwise pressure distributions at each 
spanwise station, together with calculated values of section lift coef- 
ficient, interpolated for the same spanwise station, are shown in 
figure 8. The calculations for the flaps-neutral and full -span-flaps 
configuration were made according to reference 1 and those for the 
partial-span-flaps configuration were made as described herein. For 
the experimental data, jet-boundary and stream-angle corrections have 
been applied to the angle of attack, but no corrections have been applied 
to the values of c n for the effects of the model supports and the 
boom containing the pressure tubes. Although some of the disagreement 
shown may be attributed to the assumptions involved in the calculations, 
most of it is believed to be due to experimental inaccuracies both in 
the two-dimensional data used in the calculations and in the three- 
dimensional data shown herein, inasmuch as similar disagreement is 
evident for the full-span configuration as for the partial-span configu- 
ration. In general, however, the trends indicated in figure 8 serve to 
substantiate the method used to alter two-dimensional data. 


Wing Characteristics 

Numerical integration .- In the numerical integration described in 
reference 1, the multipliers T) m , T) ms , a m , and a rnfl were determined 

by harmonic analysis. For example, it was shown that 

1 1 ifHf) - 2 f if l £ - f <«> 

x m=l m 

it nut 

where ^ sin ~ was designated as ifon. By inspection, it may be seen 
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that equation ( 19 ) could also be obtained through the application of 
the trapezoidal rule (reference 3) where f(2y/b) m sin (nut/r) is the 
function of nut/r and the interval between points is it/r. Although 
equation ( 19 ) may be expected to give good results when f(2y/b) is 
approximately elliptic, in general the application of Simpson's 
(parabolic) rule to this problem has been found to give better results 
for a wider variety of curve shapes. Since r' was originally assumed 
to be even, the interval m = 0 to r may be divided into r/2 regions, 
the area of each of which is 


it 

3r 


f 



sin 


irt 

r 



sin 


(i> D* 

r 



sin 


(i ± jh 

r 


.( 20 ) 


For the entire interval m = 0 to r, 




- = 2 




( 21 ) 


where rj m is defined and used herein as 


it (T / N m] . mit 

%-Srll- ein T 


( 22 ) 


The relations between n , ri , a , and a are the same as in 

■m ' •mo*' m ' m a 

reference lj namely. 


’nr ‘ms m 


ma 


rj = 2 t) 
ms ‘m 


( m * 1) 


(23a) 


r\ = T] 
‘ms *m 


(m = |) 


( 23 b) 


'm 


_ % 


nut 


cos 
2 r 


(24) 
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a ma - 2a m 


nut 

= ti cos — 
‘m r 


(25) 


Values of ri , 
*m' 


V , a , and a are given in table V for r = 20. 
s' nr ma ■ 


For cases for which the end of the flap or aileron is at a point 
where m is not an even integer, a special form of equation (20) is 
required which gives additional multipliers for numerical integration 
in the vicinity of the end of the flap or aileron. Such a case is illus- 
trated in figure 9 and the individual areas are 

, (26) 


Area (l) = — 




2s + 3b) - s(s + 3t) - . 

TTTtT f i + 5t f 2 


6(s + ,t)t 3 


Area (2) 


n 

r 


_t3__ 

8 s(s + t ) 


. t(t + 3 s) 
+ ' 6b ' 


+ 


t( 2 t + 3 s) 
6 (s + t) 



(27) 


Area (l) + (2) = 


"fs. 


t)(s + t) (s + t)3 ( 2 t - s)(s + t) _ 

f i + 7Zt f 2 + ZtT r 


6 st 


VC 


( 28 ) 


where s and t are expressed as fractions, of Jt/r. It can be seen 
that equation (28) is the same as expression (20) when s = t = 1. An 
example of the computation of these special multipliers is given in 
table VI for 0.5-span flaps. For the area B in the sketch accompanying 


this table, equation (28) is used with s = 1 and t = and for the 

2 ^ 

‘area C, equation (26) is used with s = — and t = 1. Values of the 


‘ms 


multipliers are given in table VII for various locations of the end 


of a flap . Values of the multipliers 




nr 


and 


a can be readily 
ma 


obtained from the values of 


^ms 


by means of equations (23) to (25). 


It should be noted that two multipliers are given at the end of the 
flap. For distributions which are continuous through the end of the 
flap, such as the distribution of c^c/b, the sum of the two multipliers 
may be used. For distributions which are discontinuous at this point, 
such as the distribution of c^.c/b, each multiplier should be used 

1 2v* 

separately with its appropriate value. For example, for -*j— = 0.5, the 
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multiplier -0.03023 should he used with the value at the outboard end 
of the flap, and the multiplier 0.03930 should be used with the value 
at the inboard end of the unflapped span. 


Wing coefficients . - The formulas derived in reference 1 are 
repeated herein for convenience. For asymmetrical distributions 



For symmetrical distributions 



For asymmetrical distributions 



(29a) 


(29b) 


(30a) 


For symmetrical distributions 



For asymmetrical distributions 



For antisymmetrical distributions 

£-1 



(30b) 


(31a) 


(31b) 


For asymmetrical distributions 


' n i 


*A 

IBo 


r-l 


m=l 


c z c 

b 



(32) 
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The values of tj^, 1 \ as ^ a m > an ^ a ma should he used as defined herein, 

and the summations should also include the special multipliers in the 
vicinity of the end of the flap or aileron. 


ILLUSTRATIVE EXAMPLES 
Symmetrical Distributions 


The method described is applied herein to a wing, the geometric 
characteristics of which are given in table VIII. The section lift 
curves are shown in figure 10 as dotted lines before being altered and 
as solid lines after being altered. The unaltered data are obtained by 
interpolation of two-dimensional data for the proper Reynolds number and 
airfoil thickness -chord ratio for each spanwise station in a manner 
similar to that of reference 1 . The altered curves’ are obtained as 
described previously and are used in the calculation of the lift dis- 
tributions in tables IX and X. These tables were designed for use with 
calculating machines capable of performing accumulative multiplication. 
The mechanics of computing are explained in the tables, but a few items 
need additional explanation. 

4 > 

The initial approximation of c z c/b is obtained in table IX. The 
computations in this table are based on the method of reference 6, the 
pertinent equations of which axe modified to suit the present purpose . 
The spanwise stations, at which values of c z c/b are computed, are 
listed in column (1). Columns (2) to ( 4 ) axe used to obtain the addi- 
tional lift distribution c Za ^ c jc for Cl = 1 according to the 

approximate relation 



Columns ( 5 ) to (7) are used to obtain faired values of the twist €5 
due to flap deflection according to the relations 




(3ta) 
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This fairing is illustrated in figure 11. Columns (8) and (9) are used 
to obtain the average twist according to the equation 



(35) 


The value of. obtained by this summation is given at the bottom of 
the table. • Columns (10) to (12) are used to obtain an single of attack 
a' defined, for the flapped span, as 


a ' = AE I \ ( g 6 " € b) + e 6 " 6 + a * (36a) 

and, for the unflapped span, as 

a ' = AE + 6 ( e & _ ' € b) + e & + a 


The values of the geometric angle of attack a' in column (ll) include 
the values of continuous twist €. Values of ij in column (13) 

are obtained for the section lift curves at the angle a'.. In t^e non- 
linear range of the section lift curves, different values of Cj ( a i) 

are found corresponding to the two values of a' at the end of the 
flaps. For the purpose of this table, the average value of c ^ ( a >) is 

used. Finally, columns (lk) and (15) are used to obtain the initial 
approximate values of cjc/b according to the relation 


( c l c \ C *al C 

W Approx ' =<“ + 2 > Cl(a ') 

which was adapted from the equation 

c l c a g Cl al C 

b ~ AE + 2 c 


s + 


AE_+2 
AE + 6 


( 6 - € ) 


(37) 


(38) 


given i n reference 6. 
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Tjje value of 8 used in table IX should correspond to the value 
of Ci obtained at the bottom of the table. Some value of 6 must 
be selected, ^however, for the computations needed to determine the 
value of C7 . It is therefore advisable to omit the computations in 
columns (10) to (15), except the bottom row, until a value of 8 has 
been found by trial and error which is consistent with the value of c 1 . 


If the lift distribution is to be obtained at a relatively high 
angle of attack, as in the example shown, the values of a' may be in 
the rapidly curving ^art of the section lift curves or even greater than 
the angle of attack for maximum section lift coefficient; For the pur- 
pose of this table, the section lift curves may be extended in the gen- 
eral direction of the nearly linear part of the curves in order to 
obtain the values of Cj( a i). This procedure is justified inasmuch as 

the subsequent operation in column (15) will reduce the value of c^. 


The approximate values of cjc/b (column (15), table IX) are used 
in column (3) of table X and the computations indicated are performed. 
Except for the computation pf a c , the computa'tions in this table are 
the same as those in reference 1 for wings without flaps. The value 
of 8 used in column (16) is obtained at the section lift coefficient 
c 1 corresponding to the value of (c^c/b)* obtained by means of the 

interpolation multipliers v m . • 


In this method of successive approximations, the value of c^c/b 
in column (22) will usually not check the initial approximate values in 
column (3)- The values to be used in subsequent approximations may be 
found by the equation 



(39) 


where A' ( c jc/b) is the increment to be added to the approximate values 
to obtain succeeding approximate values, A(cjc/b) is the difference 
between the check values and the approximate values (column (22) minus 
column (3) ), and K and Ki are constants for any particular wing. 
Equation (39) is derived in the appendix, and values of K and Ki 
for r = 20 are presented in figure 12 as functions of AE/t). These 
values compare favorably with those empirically determined in refer- 
ence 1 (Ko - 3, % - 1, K t>1 . 0 , and K = 8 to 10 j . Although the se 

values were obtained for elliptic wings, they can he used for wings of 
other plan form. The number of terms of equation (39) needed for any 
particular approximation depends upon the convergence of the approxima- 
tion; fewer terms are needed when the differences A(c2 c/b) are small or 
when positive differences nearly cancel negative differences. Values 
of Kj^ for values of i greater than 3 are small enough to be considered 
negligible. 
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Some additional explanation of the use of equation ( 39 ) may be 
desirable. For values of m from 3 to 7, equation (39) may be used 
directly if the constants ^>3 Eire neglected. For values of m 

from 8 to 10 , equation ( 39 ) may be expanded as 


*0 i ( f f ) 8 + ( K 1 + K a) 4TI * *2 A (~r) 10 


( hOa) 


££) 9 ‘ i k 3 a (T ) 6 * *a A ( ! r ) 7 + ( K i * k 3> 


(Ko * K 2 ) ♦ Kl a(^) 


V 0 '10 


(40b) 



(40c) 

for symmetrical distributions, since a(^-J = A Ir^-) . For values 

' h > 'r-i 

of m of 1 and 2, equation ( 39 ) may be expanded as 
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since 



After convergence is obtained in table X, the values of cjc/b 
(column (3)) and aq (column (l8)). are entered in table XI and the lift 
and induced-drag coefficients axe determined through use of the appro- 
priate multipliers tj^ . 


Asymmetrical Distributions 

If the angle -of -attack distribution is not symmetrical, the asym- 
metrical multipliers (3^ must be used in the nonlinear range of the 

section lift curves. Typical asymmetrical distributions are those for 
a rolling wing or for a wing with deflected ailerons. Illustrated in 
table XII is the case of a wing with flaps but without deflected ailerons, 
which is rolling at such a rate that the tip helix, angle pb/2V gener- 
ated is 0.01 radian or 0.573°- Added to the normal angle of attack 

is an increment equal to ^ ^ (in deg), which is the equivalent 

twist of the rolling wing. In order to reduce the size of the computing 
form, the table of the multipliers (3^ is arranged in the form shown. 

The values of (3^ for a positive value of 2y/b are the reverse of 
those for the corresponding negative value. Instead of reversing the 
values of in table XII, the values of c^c/b are written in 

reverse order in column (l4). Using these values with the values of 
ft mk gives the uncorrected angle 0 ^ for the stations listed at the 
bottom of the table. 

In general, the value of 8 will be different for the two sides 
of the wing,' so that the appropriate values must be used with columns (l6) 

and (17)- The values of cc c /8_ in column (l6) are those for = -0.6 

from table II, while the values of a c /5 + in column (17) are the 

negative of those for = 0.6. 

b 

Another modification must be made for asymmetrical lift distribu- 
tions since the edge-velocity factor E' should be used for the antisym- 
metrical part of the distribution ( see reference 6) . whereas the edge- 
velocity factor E has been used to alter the two-dimensional lift 
curves. From reference 6 for unswept wings. 



(42) 


E ' 
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This ..value may be taken into account in the following manner: The 

symmetrical part of the effective angle distribution is 

(a - ai) k + (a - a,i) r _ k 
2 

which is used directly with the altered lift curves. The antisymmetrical 
increment in the angle distribution is 

(a - ai ) k - (a - 

which must-be multiplied by the ratio E/E' in order for it to be used 
with the same lift curve. The effective angle is therefore 

(a - a ± ) k + (a - a 1 ) r _ k £ (a - a i ) k - (a - a i ) r _ k 

■“» = 2 + I 7 “ 2 ~ 


m 


(44) 


for table XII 


E' - E 


» (a - 04 ) - 

1 k 2E 


= a - - AXg 


(a - 04) - (a - 04) 




where 


Ai e = 


E' - E 
2E 1 


(a- - 04 ) - ( 

k 


a - “I’r.kJ 


Equation (44) is computed in column (.21) of table XII. 


Other than these modifications, the computing required 
is similar to that for table X. 


DISCUSSION 


The lift characteristics of two wings without flaps and with 
60 -percent and full-span flaps have been calculated by the method 
described herein and are presented in figure 13 together with expert- ' 
mental results from reference 7 . One wing had NACA 64-210 sections and 
was equipped with split flaps. The other wing had NACA 65-210 sections 
and was equipped with split, single slotted, and double slotted flaps. 

For the split-flap conditions, the agreement between calculated and 
experimental results is quite satisfactory, whereas the agreement is less 
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satisfactory for the single- and double -slotted-f lap conditions. Since 
the discrepancies occur for both the 60-percent-span and full-span con- 
ditions, they are probably due to differences between the two-dimensional 
and three-dimensional section characteristics rather than due to the 
method of calculating. Some of the discrepancies in maximum lift coef- 
ficient may be due to the fact that the characteristics of these wings 
were extremely sensitive to small surface irregularities. 

The stalling characteristics of these same wing-flap combinations 
are presented in figure 14 together with the calculated stall- margin 
distributions. The stall margin Ac 1 is the difference between the 
maximum section lift coefficient altered as described herein and the 
section lift coefficient at the maximum wing lift coefficient. The 
spanwise location of zero margin should correspond to the location of 
initial stall, and the margin at other spanwise locations is an indica- 
tion of the manner in which the stall spreads. In general,' the agree- 
ment between the experimental and calculated stalling characteristics 
is very good. 

The foregoing comparisons between calculated and experimental 
results were made for the same Reynolds number. If possible, such com- 
parisons should be made at the same Mach number also, unless the Mach 
number is low enough to have, a negligible effect. Even at relatively 
low values of free-stream Mach number, adverse compressibility effects 
on maximum lift coefficient have been noted (reference 8) when sonic 
velocity is reached locally on a wing. Similar effects in two- 
dimensional flow have not as yet been thoroughly investigated so that 
calculations based on available two-dimensional data must be limited to 
subcritical Mach numbers. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , November 13, 1950 
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APPENDIX 

DETERMINATION OF COEFFICIENTS USED TO 
OBTAIN SUCCEEDING APPROXIMATIONS 


In the method of successive approximations to determine the lift 
distribution, it is desirable to reach convergence with a minimum number 
of approximations. This desideration necessitates that each successive 
approximation be obtained from preceding computations in some manner. 

The manner in which these operations were performed in reference 1 and 
the coefficients used therein were determined empirically. It is shown 
hereinafter, however, that the procedure and. similar coefficients can 
be determined theoretically. 

In the following derivation, A(cjc/b) is used to designate the 
difference between the check values and the approximate values, and 
A'(c^c/b) is used to designate the increment to.be added to the approxi- 
mate values to obtain the succeeding approximate values. The section 
lift curves are assumed to be linear and 5 to be zero. The operations 
performed during the first approximation (table X) may then be repre- 
sented by the equation 





J 


If suitable increments A'(c^c/b) are chosen so that the check values 
become equal to the approximate values for the second approximation. 



The difference between these equations is 




l80r 


8k 



(A4) 


From reference 1 
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The lift -curve slope may be expressed as ' 


= 3 . — 
a o " E 90 


(A5) 


where . t| is the ratio of the actual two-dimensional lift-curve slope 
to the theoretical thin-airfoil value. For elliptic wings. 


'c \ 4 . k« 

PJk «A r 


(A 6 ) 


Therefore 


— ) X - m 

b / k ^kk - AE 


(A7) 


which, for constant q, is constant for all spanwise stations, 
stituting this value into equation (A3) yields 


Sub- 


k mf k . m k 

Equation (a 8 ) represents r /2 simultaneous equations which may be. 
represented in matrix form as 


M 

C 7 C 

A 1 -i- 

sal 

11 

C 7 c 

A -4- 


b 

qr 

b 


(A9) 


where [jlj is a matrix with all the principal diagonal elements equal 

to 1 + — and the other elements are Vik/^kk' This matrix can 

readily be put into a symmetrical form and its reciprocal obtained by 
one of the standard methods presented in reference 9 . Then 


c z c 


AE 


M 


t] 


(A10) 


bj T|r | | 

For a given value of AE/q, equation (A10) may be expressed in the form 


1-1 
2 


*(¥) + . 

N m . m±i 


(All) 
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because of the particular properties of this reciprocal matrix. For 
convenience K]_ is made equal to unity and the values of K and % 
are adjusted accordingly. The values given in table XIII and figure 12 
were obtained for various values of ' AE/q and r = 20. 
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TABLE V.- WING- COEFFICIENT MULTIPLIERS 
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TABLE VI.- CALCULATION OP AREA MULTIPLIERS UtTTH SPECIAL 
CONSIDERATION AT THE END OP A 0. 5-SPAN FLAP 



0.156 0.509 o.l*5k 0.588 0.707 0.809 0.891 0.951 0.988 1.0 
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TABLE EX.- CALCULATION OF INITIAL APPROXIMATION OF LIFT DISTRIBUTION FOR EXAMPLE WING WITH 0.6 SPAN FLAPS 
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Operations indicated in column headings for 0 < =£- < operations at bottom for 



TABLE X— CALCULATION OF LIFT DISTRIBUTION FOR EXAMPLE WING WITH 06 SPAN FLAPS 
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TABLE XIII.- COEFFICIENTS USED TO OBTAIN 
SUCCEEDING APPROXIMATIONS 


AE 
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Ko . 

K l 

*2 

k 3 
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2.384 

1.000 
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0.339 
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8.418 

3.000 

1.000 

.401 

.258 

12 

8.U1 

3.571 

1.000 

• 338 

.215 




Figure 1.- Typical lift distributions for a wing with 0.5-span flaps 









Figure Values of the lift distribution per degree of discontinuity in 
induced angle of attack as a* function of the spanwise location of the 
discontinuity 2y*/b. - 




Figure. 5.- Typical curves used to obtain factor for altering two-dimensional 

data. 







Figure 6.- Factor for altering two-dimensional data for several wings. 
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(calculated) 
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Figure 8.- Continued. 
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(d) Full-span split flaps 



Figure 9-- Definition of symbols used in numerical integration for unequal 

increments in independent variable. 
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Figure 10.- Unaltered two-dimensional and altered section lift curves for 
example wing with 0.6-span flaps. 
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<r , deg 

(c)> NACA 65-210 sections; single slotted flaps. 
Figure 13.- Continued. 
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<*L 



ac , deg 


(d) NACA 65-210 sections; double slotted flaps. 
Figure 13.- Concluded. 










Plain wing 0.6-span flaps | Pull-span flaps 



(b) MCA 65-210 sections; split flap 


















Plain wing 0.6-span flaps Pull-span flaps 



(d) NACA 65-210 sections; double slotted flaps. 






